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Abstract 
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whose asymmetry could reach 10 — 15% in certain parts of the parameter space. 
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I. INTRODUCTION 



Discrete charge-parity (CP) symmetry violation is an interesting and puzzling phe- 
nomenon within the standard model (SM). While the SM has been successful in explaining 
all of the available experimental data, it fails to provide an explanation for the breaking of 
CP. Insight into understanding its nature and structure would shed light on diverse issues 
ranging from the origin of the mass to the evolution of our universe. 

The SM presents an economical scenario for CP violation, described through only one 
weak phase in Cabibbo-Kobayashi-Maskawa (CKM) matrix. Even though the observed 
direct CP violation in K [l| and B 2] decays can be accommodated via the CKM matrix 
of the SM, this is not true for other phenomena, such as the matter-antimatter asymmetry 
present in the universe. In addition, the Standard Model fails to explain the origin of CP 
violation but merely parameterizes it. Thus the SM framework could be viewed as a low 
energy effective version of a more complete theory, and CP violation offers a motivation 
to go beyond the SM. In the models beyond the SM, the existence of new sources of CP 
violation other than the CKM phase, resolves some of the problems of the SM, but not all. 
One of the leading candidates of physics beyond the SM is supersymmetry, in particular the 
minimal supersymmetric standard model (MSSM). Supersymmetry provides a compelling 
argument for the stabilization of the Higgs sector against quadratic divergencies and allows 
unification of gauge couplings at high energies, both unexplained by the SM. In the MSSM, 
enough baryon asymmetry can be generated at the electroweak scale at any temperature by 
means of the existence of an additional Higgs doublet. However, the MSSM has difficulty 
reconciling the smallness of the electric dipole moments with expectations of scalar fermion 
masses and the size of new CP violating phases introduced by soft supersymmetry breaking 

In MSSM, there are many parameters which can in principle be complex, even after 
making all allowed rotations to get rid of unphysical phases. This raises the interesting 
possibility that CP violation, or CP asymmetries, might arise in other sectors of the model 
than the quark sector, and that would provide a spectacular signal of physics beyond the 
standard model. In particular, it is possible that there might be a close relationship between 
the Higgs sector and CP violation. The structure and properties of the Higgs sector are 
under intense scrutiny at present. Though indications at LEP for a light Higgs boson of 
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mass around 115 GeV are encouraging, they await confirmation. 

Searching for Higgs bosons is one of the major objectives of present and future high 
energy experiments. In particular, a charged Higgs boson, predicted by most models to 
have mass of the order of the weak scale, would be a definite sign for physics beyond the 
SM. The phenomenology of charged Higgs bosons is less model dependent than that of 
neutral Higgs bosons, and it is governed by the values of tan /3 and ra#± . Because charged 
Higgs couplings are proportional to fermion masses, the decays to third generation quarks 
and leptons are dominant. At hadron colliders, such as the Fermilab Tevatron and the 
CERN Large Hadron Collider (LHC), a light charged Higgs boson can be produced from 
the decay of top quark via t — > H + b, if m H ± < m t — mi,. If the charged Higgs boson is 
heavier than the top quark, there are three channels for producing charged Higgs pairs: 
pp — > H + H~, pp — > W^H^ and gb — > tH~; as well as the single charged Higgs production 
cs,cb — > H~; see {4] and references therein. In many cases complementary information 
from more than one channel will be accessible at the LHC. The LHC has a high potential 
for detecting heavy Higgs states which might be beyond the kinematic reach of the planned 
International Linear Collider (ILC). At the ILC, the main production mechanism for charged 
Higgs bosons would be e + e~ — > H + H~ , followed by one of the allowed decays H + — > tb, tv t 
or cs. The pair production cross section for charged Higgs at ILC is about 2.5 larger than 
for the neutral ones {5J. Provided that a Higgs boson couples to the Z boson, the ILC will 
observe it independently of its decay characteristics. The discovery potential is practically 
independent of tan (5 and extends up to 1.2 TeV for an integrated luminosity of 3000 fb _1 [5[. 
For specific examples of how the integrated information obtained by ILC and LHC can be 
used for Higgs detection and determination of parameters in the Higgs sector, see Q]. 

In its most general form, the MSSM predicts a plethora of new CP phases. These new 
sources of flavor and CP violation give rise to the enhancement of CP violation effects alluded 
to before, which could provide distinguishing signs for MSSM at present and future colliders. 
In this study, we concentrate on analysis of the CP asymmetry in charged Higgs decays 
H~~ — > Uidj in the framework of the MSSM. Here Uidj = cb, cs, ib, is. The corresponding 
neutral Higgs CP asymmetry in h — > didj has been discussed in Ref . [t| • The CP asymmetry 
in the main decay mode H~ — > ib and other two-body non-quark charged Higgs decays have 
been considered in jg, 9]. Recently, these authors have discussed various CP asymmetries in 
H~ — > ib by including the decay products of the top quark and subsequently the W boson 
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and showed that the decay rate CP asymmetry can go up to 25% [10J. In this work, we 
revisit the asymmetries in H~ — > ib, but also discuss the other three quark decay modes and 
compare the size of the CP asymmetry in all the channels. We show that even though in 
some part of the parameter space, the is channel has sizable CP asymmetry with respect to 
the ib, cb, or cs channel, this result has to be taken with caution since the former has very 
small branching ratios (Br's) which makes harder to observe it. We introduce and discuss 



another relevant quantity 1 ll|, {A 2 CP x Br)' 1 , for each decay mode, with Aqp the CP 
asymmetry, and Br the branching ratio for a given decay mode. This function was shown 
to be a measure of the number of required charged Higgs bosons to be produced at colliders 
for observing an asymmetry for a given channel. Based on this analysis, we conclude that 
H~ — > cb, H~ — > cs, and H~ — > ib are all optimal channels which could reveal a measurable 
CP asymmetry at the order of 10 — 15%. We also discuss the CP asymmetry induced by the 
phases of the flavor violating parameters 5lf D alone and note that their contribution is small 
and thus unlikely to account for a measurable CP asymmetry in any of the decay modes 
considered in this study. 

The outline of the paper is as follows. In the next section, we briefly review the MSSM, 
concentrating on the sources of CP violation. The decay processes are presented in Section 
III 1 1 and the numerical analysis of the decays under consideration in the Section IIVI We 
conclude in Section [V] Details of the method used for calculations are presented in the 
Appendix. 

II. THE UNCONSTRAINED MINIMAL SUPERSYMMETRIC STANDARD 
MODEL 

The superpotential W of the MSSM Lagrangian and the relevant part of the soft breaking 

T • /^squark , ■ 1 

Lagrangian L s ^ ft are respectively 

w = fiH 1 H 2 + y; ! iri:< j l; + y;/ n x (j ( p u + v;'7/ 2 gv;, (2.1) 

Cr k =-Q 4t (M| W - <r(\if.),^ - &(M*) V $ + Y^Q.H'uj + Y^Q.H 1 ^, (2.2) 



1 {Agp x Br) -1 is closely related to the total number of events TV required to establish a measurable CP 
violation for a particular mode. The exact formula is N = s 2 (Aq P x Br e) _1 . Here s is the standard 
deviation and e is the detection efficiency. 
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where H l and H 2 are the Higgs doublets with vacuum expectation values v\ and v-i re- 
spectively, Q is the 57/(2) scalar doublet, u, d are the up- and down-quark 57/(2) singlets, 
respectively, Q, it, d represent scalar quarks, Y U)d are the Yukawa couplings and i, j are gener- 
ation indices. Here A^ represent the trilinear scalar couplings. In Eq. (12.21) we are assuming 
a chiral limit of MSSM. 

We work in the unconstrained version of the MSSM and use the mass eigenstate method 
Q, where 8qUark m as 8 matrices a r e given in the super-CKM basis, and a r e diagonals 
by rotating the superfields. In this basis, the up-squark and down-squark mass matrices 
are correlated by this rotation and thus not independent. Potential new sources of flavor 
violation arise from couplings of quarks and squarks to gauginos. This method has the 
advantage that, when the off-diagonal elements in the squark mass matrices become large, 
the method is still valid, unlike perturbation-based expansions. The up(down)-squark mass 
matrices between second and third generations are taken as 
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where we assume a general flavor violation among families. In addition to the flavor depen- 
dence in Yukawa matrices, there are additional sources of flavor violation due to the soft 
mass terms and A-terms. The richness of the flavor structure depends on the assumed tex- 
tures of soft mass and A terms at the GUT scale. Assuming both the soft mass and A-terms 
to be flavor blind at GUT scale can still induce flavor violation at electroweak scale due 
to their evolution from GUT scale down to the electroweak scale using the renormalization 
group equations. The flavor violating effects become much big ger if the soft mass and/or 



74-term(s) are flavor dependent at GUT scale (see, for example, 



13j for details). 



For reasons we discuss further on, we assume fj, real and a common phase for A C)t and 
A Sy b- Note that A s has negligible effect since it is multiplied by the strange quark mass, 
tan (3 is the ratio of the vacuum expectation values of the two neutral Higgs bosons. 

From the mass matrix, Eq. (12. 3p . we also define the flavor mixing parameters as scaled 
off-diagonal flavor violating entries 



! M 2 _ \23 

(*u(d))ab= * , {A,B = L,R), (2.5) 

iw SUSY 

where Msusy is the common scale for the parameters M~ and M^,^ g . As mentioned 
before we allow S^r D ^s to be complex (to have CP violating phases). 

We do not repeat listing the chargino and neutralino sectors of the MSSM here, as we 
don't assume a new non-zero CP phase in either sector; for the details see Ref. jl^j. As the 
gluino contribution is dominant for the charged Higgs decays, we give the relevant up-type 
quark-squark-gluino interaction g: 

3 

C U u 9 = J2^9s T r st [fi« (IV)* Pl 9 r fi* - u\ (IV) (<+8)a P R g r u\ + H.c] , (2.6) 

i=l 

where T r are the SU(3) C generators, Pl,r = (1 T 7s)/2, i = 1, 2, 3 is the generation index, 
a = 1, . . . , 6 is the scalar quark index, and s, t are color indices. There is a similar interaction 
for the down case. 

Considering the number of complex parameters in MSSM, further assumptions are needed 
for simplicity and predictability. In the gaugino sector, two of three gaugino masses could 
be complex, unless a degenerate spectrum at grand unified theory (GUT) scale is assumed. 
We assume all three gaugino mass terms real at any scale. The higgsino mass parameter \i is 
in principle complex. However, the phase of a is strongly constrained by the electron dipole 



moment (EDM) measurements of neutron 15| and cannot exceed values of the order of 0.0 



0.001. For the effects of CP violating SUSY phases on other EDMs and systems, see [16 ] . 
We simply neglect this phase and consider /x real. In the squark sector, there are trilinear 
soft couplings of quarks, A u ^, which are complex. In addition to these, the misalignment 
between quarks and squarks arising through the diagonalization of their respective mass 
matrices leads to new sources of flavor violation, with parameters denoted as 5u,d- These 
are also generally complex, and could have large imaginary parts. We will assume non-zero 
flavor violation only between the second and third generations, because the ones involving 
the first generation are required to be small, based on experimental constraints in K and 



D physics {13I . \v\ . The kaon mass splitting parameter AMk due to K — K mixing, and 
the parameters e and e' put severe constraints on the real and imaginary parts of the flavor 
violating parameters 5^ 2 , respectively. Somewhat weaker bounds can also be obtained from 
K physics for the up sector between the first and second generations as well but stringent 
bounds on 5]f can be obtained by using the experimental bound on the D — D mixing 
parameter AM p. In a similar fashion, the flavor violating parameters 5^ between the first 
and third generations can be restricted with the use of low energy B physics measurements 
(AM Bd ,B —> X s j, Sb^k s )- The common feature of these constraints is that the upper 
bounds on the flavor parameters involving the first generation have to be less than 0.1 
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17| . There are however no similar limits on the mixings between the second and third 



generations so that we will keep them arbitrary. 

We assume that there are non-zero phases from the trilinear couplings A u ^ and intergen- 
erational flavor mixing parameters 5^ D . For further simplification, we assume a common 
phase ArgfAJ = Arg[A d ] = ArgL4]. 

Thus, the supersymmetric sources of CP violation of interest in charged Higgs decays 
come from the soft broken terms M?, M?. and the trilinear scalar coupling A u , and are 
introduced through the matrix Ty. In the following section, we analyze their effects on the 
calculation of the CP asymmetry. 



III. CHARGED HIGGS DECAYS H~ -» u { dj 



In this section we discuss the CP asymmetry in the charged Higgs boson decays H — > 
Uidj, which is defined as 2 

CP T(H- -> Uidj) + T(H+ -> Uidj) ' 1 ' ' 

where V is the partial decay width of the decay mode considered. These decays are tree level 
processes and one could calculate the branching ratios, decays widths, etc. with tree-level 
approximation. 

However, it is known that [3] the CP-odd observable Acp requires a nontrivial phase 
from Feynman diagrams (called absorptive or strong phase), in addition to the weak phase 



This type of CP asymmetry is sometimes called partial rate asymmetry. For other types see 10|, [18 1 
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FIG. 1: The tree and relevant one- loop diagrams contributing CP asymmetry for the decays 
H~ — > Uidj. 

mentioned in the previous section. This way, the imaginary part of the amplitude is non-zero, 
Zm (Amplitude) ^ 0. One way of introducing such a phase is through one-loop Feynman 
diagrams, where some of the intermediate particles go on-shell. Then, the numerator of 
Eq. (13. II) will be proportional to the interference term between tree level and one-loop 
contributions. The tree level and relevant one-loop contributions to the decays H~ — > Uidj 
are shown in FigdJ 

There are many more one-loop contributions to the decays, but based on various kinemat- 
ical considerations only six of them are relevant to Ac p (four vertex type and two self-energy 
type diagrams) 3 . We consider cuts through chargino-neutralino internal lines, or through up 
squark-down squark lines. We will call them internal- cut states. They contribute to the CP 
asymmetry when ran ^ m^o + and m# ^ + m d - are satisfied, so that the chargino 
and the neutralino can go on-shell. This is the necessary condition to induce an absorptive 
part into the amplitude. The generic self-energy and vertex type diagrams with cuts are 
shown in Fig. [2] where the internal-cut states represent two possibilities in each case. Here 

3 There exist additionally some SM vertex contributions with W boson, CP-even Higgs bosons, and quarks 
in the loop. We don't get CP asymmetry contributions from such diagrams since we take fx real. For the 
case with complex u, see Q. 
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we don't count the nipping cases (k\ <-» A^), which are only relevant to the vertex type 
diagrams, in which case the third intermediate state is different. For the vertex-type dia- 
grams, Fig. &[b), there are two other possible cuts, through k\ — k 3 or k 2 — k 3 . The cut with 
k\ — &3 is not kinematically allowed since it requires ^ m d a + m Xn > md j ^ m d a + m x + ' 
or m^. ^ mj + m g . None of these conditions is possible. 



Pi 
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(a) 



(b) 



FIG. 2: The unitarity cuts for self energy (a) and vertex type (b) generic diagrams. 



For the case with ki — k 3 cut, based on the current experimental bounds on the mass of 
the lightest neutralino and the lightest up-type squarks, there could be some contributions 
to the decay mode H~ — > idj from the diagrams (2) and (4) of Fig. [JJ in a very narrow 
kinematical range (m„ + m^o < m t ). We simply neglect such contributions. There exists yet 
another way to produce the necessary absorptive cut, by taking the invariant mass squared 



for the final states, {j>2 + P3) 2 , to be greater than (m^o 



which results in some cuts 



19j . in the analysis of the three-body 



in the phase space. This method was persued in 
semileptonic top quark decays. 

One can show that the numerator of Eq. (13.11) is proportional to the imaginary part of the 
amplitude from loop diagrams 4 arising from tree level-loop interference terms. As usual, we 
neglect possible loop-loop interference effects which are much smaller than tree-loop terms. 
In Appendix, instead of giving rather lengthy analytical results for this imaginary part, we 
outline the method used to do the calculation using the cuts in Fig. [2j We present it in a 



generic way, valid for the four decay modes of the charged Higgs boson, H~ 



Uidj. 



For example see 18( for the details 
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IV. NUMERICAL ANALYSIS 



In this section, we present comparatively our numerical results for the decays if" — > 
iiidj (ui = c,t, dj = s,b) 5 . In the calculation of the CP asymmetry from the formula in 
Eq ( 13. ip . we use the tree level values instead of one loop when evaluating the sum of the 
partial decay widths (T) in the denominator, which simplifies the numerical calculations. 
This approximation is not valid for the numerator where one loop contributions are needed 
to extract the CP phases. 




100 200 300 400 500 600 700 800 900 1000 100 200 300 400 500 600 700 800 900 1000 



m H -(GeV) m H -(GeV) 

FIG. 3: The partial decay widths and branching ratios of the charged Higgs boson decays H~ — > 
Uidj. For the supersymmetric channels; M\ = 95 GeV, M2 = 200 GeV, [i = 250 GeV, Msusy = 
500 GeV, and A = 400 GeV. 

5 The CP asymmetry in one of these decays, H~ — > ib has been discussed in Q under a different parameter 
set and assumptions. A good agreement is obtained once we switch our parameter values to their set. 
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In Fig. [3] we show the partial decay widths of the channels H~ — > cb, cs, ib, is at tree 
level, together with their branching fractions (Br) for small and intermediate tan (3 values. 
We represent the other tree level decay channels as u the rest 1 and we use the FeynHiggs 



program |25|J to calculate the partial decay widths of these channels. Basically, depending 
on the charged Higgs mass, "the rest 1 includes the lepton channels H~ — > ev e , [iv^, tv t , 
the neutralino-chargino channels H + — > xtxp * = 1, 2, j = 1...4, the Higgs-vector boson 
channels —>■ h°W, H°W, A°W, and the sfermion channels H ± — > fcfj i,j = 1...3. 
Inclusion of these channels is important as it affects the number of charged Higgs bosons 
required to observe the asymmetry, as discussed later in this section. 

In the region m H + ^ m t , H~ — > ib is the dominant mode. However, in the region 
m H + ^ 140 GeV, the leptonic decay — > tv t becomes the main decay channel, and 
in between these regions, including threshold effects, the below-threshold three body decay 
if ± — > h W* has a branching ratio comparable to, or even dominating over other channels, 
its exact value depending on the mixing in the Higgs sector. Inclusion of threshold effects 
also opens other three body channels, like H ± — > A W* and H _ — > bi* that have sizable 
branching ratios in the intermediate mass range. We do not include the threshold effects, but 



sec 



26( | for details. We note that for heavy charged Higgs, the neutralino-chargino channels 



are comparable with the ib channel. However, it is important to observe that among these 
decays, H~ — > cs and H~ — > cb have non-negligible branching ratios. This is important for 
our analysis, since observability of CP asymmetry in a specific channel requires not only a 
sizable asymmetry but also an experimentally viable branching ratio. We must comment 
on the strange quark mass dependence: as seen from Fig. [3l while the partial decay widths 
for ib and is channels are suppressed as tan/3 gets larger, the opposite is true for the cb 
channel. This is due to the fact that, evaluating the Feynman diagrams in Fig. [H the scalar 
quark couplings to the charged Higgs are proportional to (m u cot /S + m^tan (3), as seen from 
Eq. (12. 3p . There is no supression in the cs curve since for intermediate tan/3 values the term 
proportional to the strange quark mass is comparable to the one with charm quark mass 
and becomes dominant for larger tan (3 values. So, one must keep the strange quark mass 
non-zero, at least for cs decay. Of course, its effect in the is case is negligible. 

For the numerical analysis, we fix some of parameters of the model globally because the 
CP asymmetry is not very sensitive to their variations. As mentioned before, we introduce 
a common phase and magnitude for the trilinear couplings A u and A^ as Arg[A] and A. 
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The Higgs parameter /x is taken real. We allow the gluino mass m g to be light. We also 
use the parameterization for squark mass matrices where a common scale Msusy is chosen 
for the soft breaking parameters M~ ^ and M^.^ ? . The flavor violating parameters (5's ) 
are set to zero everywhere except for the case where we test the sensitivity of Aqp to these 
parameters. There exist two other free parameters in the Higgs sector of the MSSM, taken 
in a popular framework to be tan (3, and one of the Higgs boson masses, often taken as the 
CP-odd Higgs mass uia- Of course we could equally well assume any of the others as the 
free Higgs mass parameter. As it is more convenient for the present analysis, we choose 
m>H+ as the free Higgs mass here. Unless otherwise stated, we fix the following parameters 
globally 6 

M SUSY = 500 GeV, n = 250 GeV, tan/3 = 10, M 2 = 200 GeV, m g = 250 GeV. (4.1) 

In Fig. H] we illustrate the dependence of the CP asymmetry Acp on the phase Arg[A] 
for each decay for a variety of A values. As expected, the behaviour is periodic and the 
sensitivity to both Arg[A] and A is quite significant. It seems that each decay mode can have 
an asymmetry as large as 10 — 15%. It is also seen that the asymmetry produced in H~~ — > cs 
and H — > is decay modes could be comparable, or even larger than H — > cb or H~ — > ib. 
Of course this is not very unusual, since it is possible to get large asymmetries for decays 
with smaller branching ratios [3]. Therefore, determining the optimum channel among these 
decays requires further analysis, but qualitatively one can predict that H~ — > cb and H~ — > 
ib have similar asymmetries. The cb mode has non negligible branching ratios due to the 
fact that CKM supression is compensated by the large mj mass appearing in the couplings. 
In the parameter space that we explore, the diagrams with chargino-neutralino/neutralino 
in the loop give negligible results compared to the gluino loop diagrams. So, the main 
contribution comes from gluino vertex diagram (diagram 5 of Fig. [T]), while the self energy 
diagram (diagram 7 of Fig. HD is also important. We keep these contributions, but include 
and check the others wherever relevant. 

Additionally, we want to test whether it is possible to account a sizable asymmetry in 
charged Higgs decays by introducing complex flavor mixing parameters between second and 
third generation quarks and keeping all the other parameters real. This is a secondary 

6 The gluino mass is consistent with TEVATRON limits Q 



12 




-4 -3 -2 -1 1 2 3 4 -4 -3 -2 -1 1 2 3 4 

Arg[A](rad) Arg[A](rad) 

FIG. 4: The CP asymmetry for the charged Higgs decays H~ — > iljdj as a function of the phase 
ArgL4] for the various values of A. 

effect since, unlike the case with no- zero phase Arg[A], getting the absorptive phase through 
such parameters requires not only a chirality flip in the squark propagators but also a mass 
insertion for flavor changing as well (this is a CP breaking and flavor violating effect). 
In Fig. Owe show the asymmetry Aqp as a function of the phase Arg[(5^ 3 )ii] for various 
charged Higgs masses. The small graphs inside each graph represent the positive asymmetry 
Acp in the logarithmic scale, so that one can distinguish the curves with different charged 
Higgs masses. With the exception of the H~ — > cb decay mode, which can have as large as 
0.3% CP asymmetry for = 300 GeV, the other decays yield negligible asymmetries. 
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FIG. 5: The same as Fig. [5] but as a function of the phase of the flavor violating parameter 



Asg[{5fj)Li\ for various rriH+ values. We take A = 1200 GeV and Arg[A] = 0. The absolute value 



of {5]j)ll is set to 0.5. The small graphs inside each graph represent the positive asymmetry Aqp 
in the logarithmic scale. 

The absolute value of (S^)ll is set to 0.5, but Aqp is not very sensitive to the value of this 
parameter. We also checked the asymmetry induced by the phases of the other S parameters 
in both the up and down sectors, but they all give negligible contributions. Therefore, it is 
unlikely that a measurable CP asymmetry in charged Higgs decays H — > Uidj is generated 
only by absorptive phase from (S^)ll- 

In the next two figures, Figs. [6] and [7J we plot the CP asymmetry Aqp, for various Arg[A] 



14 



-0.04 



0.00 



X 

g -0.08 
< 

-0.1 



-0.12 



-0.14 



Arg[A] = tt/2 
Art/fA] = 7r/4 
Arg[A] = tt/8 ■ 



200 400 600 800 1000 1200 1400 1600 
A(GeV) 



< 




200 400 600 800 1000 1200 1400 1600 
A(GeV) 



I 

X 



< 



Arg[A] 


= tt/2 


Arg[A] 


= tt/4 - 


Arg[A] 


= tt/8 



200 400 600 800 1000 1200 1400 1600 
A(GeV) 



0.16 



0.14 



0.12 



0.08 




0,06 



0.04 



0.02 



200 400 600 800 1000 1200 1400 1600 
A(GeV) 



FIG. 6: The CP asymmetry for the charged Higgs decays H — > Uidj as a function of A for the 
various values of the phase ArgL4]. 



values, as a function of A and the gluino mass nig, respectively. Acp in both of the decays 
H~ — > cb and H~ — > cs in Fig. O can be as large as 14% at Arg[A] = n/2 for large A values, 
but H~ — > ib remains slightly smaller. The CP asymmetry in H~ — > is becomes even bigger, 
but this has to be taken with some care. More than 10% asymmetry could be achieved for 
each decay mode for A = 1200 GeV and if the gluino is very light (~ 200 GeV), as shown 
in Fig. [71 We note that there exists similar dependence on the charged Higgs mass. The 
CP asymmetry also changes between ±15% as we vary tan/3. The dependence is significant 
especially for light charged Higgs masses and small tan j3 values. 
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In order to test the viability of these channels for the search of CP asymmetry, we scanned 
the parameter space spanned by the most sensitive decay channels and compared the number 
of events (which can taken as the required number of charged Higgs bosons) as function of 
(Aq P x Br) -1 . In Fig. [8], we show the scatter plots for the charged Higgs decays H~ — > tiidj 
in the (Aq P x Br)^ 1 — Acp plane. These events are obtained by running the sensitive 
parameters randomly in the parameter ranges A £ (0, 1400) GeV, Arg[A] £ (— ir, 7r), rrig £ 
(200, 1000) GeV, m H + £ (200, 1000) GeV, and tan/5 £ (1,50). Note that the x-axes are in 
logarithmic scale. The branching ratios are calculated by evaluating all the other tree level 
charged Higgs decay widths with the FeynHiggs program. Again, the small graphs show 
the positive part of the asymmetry distribution in the logarithmic scale. From these scatter 



1(3 




10° 10 2 10 4 10° 10 8 10 1D 10 12 10" 10 4 10 (i 10 s 10 10 10 12 10 14 10 16 10 18 10 20 

(AjpxBr)' 1 (AjpxBr)- 1 

FIG. 8: The scatter plot for the charged Higgs decays H~ — > Uidj in the (A^p x Br)~ l — Aqp 
plane. These are obtained by scanning the sensitive parameters A S (0, 1400) GeV, ArgL4] E 
(-7r,7r), m-g e (200, 1000) GeV, m H + G (200, 1000) GeV, and tan/3 G (1,50). The x-axes are in 
the logarithmic scale. The small graphs inside each graph represent the positive asymmetry Aqp 
in the logarithmic scale. 

plots, one can make a few observations. 

Each decay channel has a possibility to develop an asymmetry bigger than 10%. In the 
case of H~ — > cs and H~ — > is, there are quite number of events which yield asymmetries 
around 15%, which is as large an asymmetry as the other two channels can attain. But the 
number of charged Higgs bosons required to observe such an asymmetry is around 10 5 — 10 6 
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FIG. 9: The same as Fig. [8] but for non-zero Arg[(Sf^)LL] phase. A is set to 1200 GeV and taken 
real. The x-axes are in the logarithmic scale. The small graphs inside each graph represent the 
positive asymmetry Aqp in the logarithmic scale. 



for is decay, which is 2 to 4 orders of magnitude more than the number needed to make 
such a measurement in the H~ — > cb or H~ — > ib channel, respectively. On the other hand, 
H~~ — > cs yields a comparable asymmetry distribution with respect to H~ — > cb, and also 
the required number of charged Higgs is similar for some part of the parameter space which 
maximizes the asymmetry. Clearly, the parameter set for the maximal scenario is different 
for each decay modes. Such a similarity between cs and cb channels is not surprising since 
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the cs decay mode can have a branching ratio comparable to, or larger than, that of cb decay 
mode, depending on the value of tan/3. In the cs channel, since the strange quark is light, 
it is difficult to distinguish it from cd channel. 

So, we can conclude that, depending mainly on the tan/3 value, the asymmetry in both 
cb and cs can be competitive with the one from ib if there are enough charged Higgs bosons 
produced at the colliders. But one needs at least 10 2 times more statistics for cb and cs 
cases. From Fig. [HI the distribution of the events for ib modes, unlike the other channels, 
are not scattered much since in almost the entire parameter range considered in the scan, 
its branching ratio remains constant. 

In the last figure, Fig. [HI we perform the same scan of the parameter space as before, but 
with a non-zero phase Arg^ 2 /)^] in the up-type squark mixing matrix, while switching 
the phase of A off. So, instead of running the trilinear coupling A and its phase, we run 
(^"u)ll and its phase in the range (0,0.5) and (— 7r,7r), respectively. It seems that the 
largest asymmetry comes from the H~ — » cb process which could be at most as large as 
0.6 — 0.8%, but at least 10 6 number of Higgs produced are required to probe such asymmetry. 
Additionally, the asymmetry is still significantly smaller than the case where the absorptive 
phase originated from the trilinear couplings. For the other channels, not only are the 
asymmetries very small, but as a consequence the number of charged Higgs required to 
observe them is enormous. 



V. SUMMARY AND CONCLUSION 

In this study, we analyzed the possibility of obtaining measurable signals for CP asym- 
metry in the charged Higgs decays H~ — > Uidj for the second and third generation quarks 
in the MSSM. Above the top quark threshold, charged Higgs bosons decay mainly to ib, but 
decays to cb, cs, and is are also relevant. The CP asymmetry of the main decay mode ib 



and other non-quark charged Higgs decays have been considered previously {si, [I], 10]. Here 
we analyzed, discussed and compared all significant quark channels to pinpoint which one 
is more likely to produce a visible asymmetry at the colliders. A non-zero CP asymmetry 
requires an absorptive phase, for which we considered possible interference terms between 
tree-level and one-loop diagrams. Then we calculated the imaginary parts of one-loop scalar 
two- and three-point diagrams by first deriving first generic analytical formulas for the dis- 
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continuity in such diagrams with the help of Cutkosky rules. We then presented the results 
obtained for the vectorial and tensorial types of integrals. Within the MSSM framework, 
we investigated the effect on the asymmetry of two relevant absorptive phases; the common 
phase of the trilinear couplings A, and the phase from flavor the violating parameter (5^ 3 )ll- 
We analyzed their effects separately and without considering any interference effects. 

Consideration of non-zero values for Arg[A] predicted asymmetries around 10 — 15% for 
each decay mode. However, by including the requirement that {A 2 CP x Br)' 1 is proportional 
to the number of required Higgs bosons, our analysis indicates that only H~ — > ib and 
H~ — > cb, cs decay processes would be likely to induce a measurable CP asymmetry, with 
the requirement that at least 10 2 — 10 4 charged Higgs bosons be produced at the colliders. 
Unlike the phase Arg[A], the phase (5 2 / 3 )ll cannot account a sizable CP asymmetry Aqp- 
In this case, only H~ — > cb can get an asymmetry around 0.6 — 0.8%, with the requirement 
that 10 6 charged Higgs bosons must be produced at the colliders, in order to translate into 
a measurable asymmetry. The fact that the arguments of A and (5 2 / 3 )ll give rise to CP 
asymmetries of different orders of magnitude justifies a posteriori neglecting interference 
effects. 
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APPENDIX: THE METHOD - CALCULATION OF THE IMAGINARY PART 

We calculate the absorptive part of the loop diagrams by applying the Cutkosky rules. 
In general, in the loop integration, we end up with a numerator with scalar, vectorial, or 
tensorial structures (and their pseudo counterparts), depending on the types of particles 
running in the loop. We first consider the scalar case and then outline the method for the 
vectorial and tensorial cases and present the results. We discuss self energy and vertex cases 
separately. 
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a. Two-Point Function: 



In Fig. 121- (a), if we assume that particles running in the loop are scalar, we get a B type 



scalar two-point Passarino-Veltman function 



201 ] . We use the convention that all external 



momenta are incoming and the diagram (a) of Fig. [21 can be expressed as 

where mi and rri2 are the masses of the particles in the loop. By applying the Cutkosky 



rules 2l|| for Fig. [21- (a) we have, for the discontinuity across this cut, 

AB (pl,ml,ml) = J d 4 k 2vr 0(fc°) 5{k\ - m?) 2vr 0(p$ - fc°) <S((p a - fcx) 2 - m^) 

. >/A(p?,m?,T7^) , 2 2 x 

= 27TZ - 2 ® [Pi ~ ( m l + m 2j J , (A. 2) 

Pi 

where k\ and p\ are the energies of the corresponding particles. ABq represent the disconti- 
nuity of B° and is related to the imaginary part of the diagram up to a factor of 2i (see 22 ] 
for example). X(x,y,z) is the usual Kallen function defined as X(x, y, z) = (x — y — z) 2 —4:yz. 
The argument of the Heaviside function in the final result ensures that both internal-cut 
states, particles 1 and and 2, go on-shell in the loop. 

Now let's assume that we have a vectorial term in the numerator. If we denote this 
integral as then the discontinuity can be written 

&B»{p\,m\,rnl) = AB (p 2 v mj, m\) ® h** 



(A.3) 

, fi = i 

where means that should be considered as a part of the integrand of ABq term. 

In a similar manner, one can calculate a tensorial type of two-point integral. Since we will 
need terms only up to second rank tensors in our calculation, AB^ U is enough to calculate 
for the discontinuity in such cases, and the result can be expressed as 

AB^(pl,ml,ml) = AB (pl, ml, mf) ® (k ft k v ) 

= K^AB (pl,ml,ml), (A.4) 



where 



K - = Diag( K °,^, I), K o = M±2^!)! , , = y , (A.5) 
3 3 3 Ap( Api 

21 



Then the imaginary part of the self energy diagrams shown in the last two diagrams of Fig. [T] 
becomes the sum of the above terms 

Xm(Self Energy) = 1 ^ (Yj AB + 3f A£? M + AB^) , (A.6) 

i 

where Yi's include all other contributions arising from the Feynman rules and the index I 
runs over loop diagrams. 



b. Three-Point Function: 

The evaluation is similar to the two-point function case but the calculation is more cum- 
bersome. We first give the result for the discontinuity in a three-point scalar integral, known 
as Co(pi iP2iP\i fn\, m|, rn\) 7 . Here 1711,2,3 are the internal masses. Again using Cutkosky 
rules, we have 8 for the discontinuity 

AC = - /"^Tre^MA;?-^ 



7r 2 J (p s + ki) 2 - m\ 

l°g(^fl e(p?-(m 1 + m 2 ) 2 ), (A.7) 



VHpIpIpD \ a + f 3 



where 



a = pj (pi + 2ml - (P2 + P3 + m i + m D) ~ ( m i ~ m l)(pl ~ pI) > 

(3 = yJxi^MM) KpIpIpD- (a.8) 

The vectorial type integrals can be calculated as follows. The extra term k± can be con- 
verted into the external momenta p^s and their derivatives with the help of the propagator 



in the denominator 



23] . For example, 



2 = \ N ((P3 + hf - ml)} ~ T—A 2- (A.9) 

For convenience, we define (p 3 + &1) 2 — m\ = D, f = 1/D, fx = log(D), and / 2 = 
D (log(-D) — 1). This way we can express the discontinuities in both the vectorial and the 



7 For simplicity, we suppress the argument of C q in the rest of the paper 



Our result is consistent with the one given in 22 1 
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2 nd rank tensorial type integrals in the following form 



Ad 



Ad = AC (Jo -> /i) , 
AC 2 = AC (/o -> / 2 ) • 



(A.10) 



Then, computing AC\ and AC2 integrals in a straightforward manner, we find 
Ad = AC + ^ V / A(^?X) log O (p? - ( mi + m 2 ) 2 ) , 

e{p\-{m l +m 2 ) 2 ) (A.ll) 



Ad 



AC n - 



Pi 



A(p?,m?,m|) log 



2 e 3/2 



Here e is the Napier's constant. Next, we take the derivatives of ACi and Ad with respect 
to ^2,3 and plug them into Eq. ( 1A.10I) . After some lengthy algebra, we get 



Ad 



2iri ®(p\ — (mi + m 2 ) 2 ) 

VKpIpIpI) 



. a + (3 
l0g, a^' 



(u - au) log ( + 2/3^ P * 



a — P 



2 pj 



- 2 " 9(P; ' (m ' +m2)2) ( ^AC S + f# AC 22 + W + «) AC 23 + AC33 



AC n 



a (3v 
Apj 



,1 u (a - P)v\ . ( a + (3\ 



AC 
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(a 2 - /3 2 ) 
+2a/3 



2 ' 8p 2 

1 



a — /3 



it \ ( a — /3 

2 1 lo s 



2p? 



2pf e 



HpIpIpD 
1 



2(u - aw)' 



log ( ) + 8/3v(w - av) 



AC 23 ={3v + — — - log £ 

2 \a — p 



AG 



33 



/'a + /3 

-Pi log 



i i i 



(A.12) 
is the 



a — (3 

where u = p\ + m 2 — m 2 , and t> = (p 2 — p\ + p 2 ) / A(p 2 , p 2 ., p 2 ) and 
metric tensor. 

It remains only to determine the coefficients, like the ones defined in Eq. (IA.6j) . by 
comparing the actual matrix elements for the diagrams in Fig. [TJ The diagrams in Fig. [T] 
and the corresponding matrix elements are generated with the software FeynArts 24( and 



23 



then we do the rest of the calculation with our own code. 
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